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We elaborate further on the evolution properties of cosmological fluctuations through a bounce. 
We show this evolution to be describable either by "transmission" and "reflection" coefficients or 
by an effective unitary S— matrix. We also show that they behave in a time reversal invariant way. 
Therefore, earlier results are now interpreted in a different perspective and put on a firmer basis. 
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I. INTRODUCTION 

A (Jordan or Einstein frame-) bounce [l| is predicted 
to have taken place in some string - or quantum gravity 
- motivated early universe models, as e.g. the pre big 
bang and ekpyrotic cases @, Q. Those models are 
often presented as plausible challengers to the otherwise 
well-confirmed inflation paradigm, and must therefore be 
able to provide an almost scale-invariant perturbation 
spectrum. 

The difficulty resides in the propagation of the spec- 
trum, usually produced in the contracting phase, through 
the bounce, as the latter often involves high curvatures 
and thus knowledge of at least some non-linear terms in 
the underlying theory. Even assuming that the pertur- 
bations can in fact somehow propagate through, i.e. sup- 
posing the first order to remain meaningful all along, a 
point quite debatable in itself (see e.g. Lyth in Refs. Q), 
it is necessary to impose a specific model to examine the 
validity of the usual assumptions that are made, either 
explicitly or implicitly, in the abovementionned frame- 
works. 

This has been done in Refs. in which we used pure 
general relativity, a Friedmann-Robertson- Walker met- 
ric with positively curved spatial sections and a scalar 
field ip to provide the matter content and lead the dy- 
namics of the scale factor a(r/) as a function of confor- 
mal time 77. We have since received inquiries about some 
seemingly unexpected features of the results, such as an 
apparently vanishing-determinant transition matrix, and 
we were questioned about the unitarity of the evolution. 

In this brief report we clarify these issues by showing 
that, since the evolution equation can be understood as 
a time-independent Schrodinger equation, its solutions 
can be interpreted as transmission and reflection coeffi- 
cients, which ensures the conservation of probability, or, 
stated differentely, the unitarity of the S— matrix that 
can be built out of the transition. Note however that 
these results stem from the mathematical properties of 
the Schrodinger equation, although the actual interpreta- 
tion of the underlying classical perturbed Einstein equa- 



tions does involve neither any physical wavefunction nor 
any probabilistic meaning. In particular, this should not 
be confused with the effective quantization of the physical 
modes, i.e. the second quantization of the semi-classical 
effective field theory which underpins inflationary struc- 
ture formation scenarios. 

In addition, we show that the transition matrix reflects 
the underlying invariance under time reversal. Finally, 
we also demonstrate that the vanishing of the determi- 
nant signaled before was to be expected and turns out 
in fact to provide a necessary consistency check of the 
calculations performed in Ref. 5j. 



II. THE BOUNCE SEEN AS A POTENTIAL OF 
A SCHRODINGER EQUATION 

We want to describe the evolution of cosmological per- 
turbations in a universe with positively curved spatial 
sections through a bounce described by general relativ- 
ity and a scalar field 0. The scale factor a depends on 
the conformal time 77. One can construct a perturbation 
function u oc $, proportional to the Bardeen potential 
$ through a function of the background^ which we do 
not need to specify here, and satisfying 



v u (v)] 



0, 



(1) 



for a wavenumber k = y/ n(n + 2) on the 3-sphere. 
By setting k 2 — » E, u — » tf> and r\ — > \/2mx/h, 
Eq. is of the form of a time-independent Schrodinger 
equation -(ft 2 /2m)d 2 ip/dx 2 + [V(x) - E] tp = in 
the variable r\ for the "wavefunction" u in the po- 
tential V u (x). Given this form, the current J — 
— (ift/2m) — vbVip*) — > — i (u*vf — uu*') is con- 

served V • J = J' = 0. 

In the case of a symmetric bounce, the potential V u (rj) 
of Eq. is nonvanishing in a finite and symmetric region 
around 77 = 0, i.e. for — ry z < 77 < x z , which we shall 
call region II, surrounded by regions I (for 77 before the 
bounce) and III (for 77 after the bounce). In both regions 
I and III, the potential being negligible, we have 
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FIG. 1: The analogy and differences between the wavefunction across a potential in ordinary quantum mechanics and primordial 
perturbations through a bounce. In both cases, the total current is conserved. 



leading to the currents 

^1,111 = 



A 



(3) 



on both sides of the bounce. Apart from the overall ir- 
relevant normalization, this is understood, in ordinary 
quantum mechanics, as the difference between current 
probabilities of wavefunctions traveling respectively in 
the — r\ (coefficient A) and +77 (coefficient B) directions. 
The transition matrix relating ( A lu , B ni ) to (A I , B l ) was 
obtained in Refs. H, @ by expanding the solution in 
the region II into an even f n (—rj) = /„(??) and an odd 
0( — v) = — 0( — functions, which can always be done 
in the case of a symmetric potential. 

Setting f u = / n ( % ), f u = d/ n (ry z )/dry, = 0fa), 
and 0' = d&(r] z )/dr], we found 



X Y 
Y* X* 



thereby defining the transition matrix T u , where 



X 



lkx ° {f u - ikf u ) (& - ik&) 



Y = - 



k (/ n0 /-0/^ 

i (/; T 0'+ fc 2 0/ n ) 



(/ n 0'-0/ n ) 



(4) 

(5) 
(6) 



both terms being calculated at the points of matching. 

From Eq. Q, one sees that the determinant of T u is 
unity, namely |X| 2 — \Y\ 2 = 1, as may be checked by 
direct calculation. This property is related to the fact 
that the propagation of perturbations through a bounce 
is similar to that of a wavefunction ip(x) through a po- 
tential V(x), as illustrated on the right of Fig. ^ If 
the incoming wavefunction has amplitude B Y , it is re- 
flected with amplitude A l: leading to a reflection coeffi- 
cient R = \A l /B l \ 2 , and the amount B IU goes through, 
so that the transmission coefficient is T = {B^/BJ 2 (as- 
suming the same wavenumber on both sides). For the 
bounce case shown on the left, the analogy makes sense 
provided one reads the figure backwards in time. Then, 



the "transmitted" wave, with amplitude A 1 , in fact the 
original perturbation, travels again in the same direction 
as the wave hitting the potential with amplitude B 
leading to a transmission coefficient T B = \A 1 /A 
while the reflection coefficient reads i? B = |-B m /.<4 m | 
Then, the current conservation, J' = 0, implies that 



in 5 
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Rt> + XL 



1, 



(7) 



as expected. Notice that this is just a rewriting of the 
matching condition |B ln | 2 - \A 1U \ 2 = |S,| 2 - |AJ 2 , itself 
a consequence of the fact that the transition matrix is 
of unit determinant: this merely expresses the conserva- 
tion of probability in ordinary quantum mechanics. The 
result obtained in Ref. is therefore in perfect agree- 
ment with the conservation of the effective probability 
current associated with the otherwise classical Einstein 
equations. 



III. BOUNCE S MATRIX AND UNITARITY 

Another way to express conservation of probability is 
through unitarity. Therefore, one must consider a differ- 
ent point of view, namely that of the S— matrix, to build 
the unitary transformation stemming from T u . This is 
what we do in this section, by means of a comparison 
with standard collision processes in ordinary quantum 
mechanics. Obviously, there is no actual collision in- 
volved, and the fact that one can construct such a uni- 
tarity S— matrix merely reflects the mathematical con- 
sistency of Refs. 

In three dimensional quantum mechanics, a plane wave 
propagating in an incident direction n = k/|k| and diffus- 
ing in a potential located at O is observed at the point 
M such that OM = rn' leads to the collisional wave 
function ^^n' given by 



4. _ / lP ikr 
Jn,n' c j 



(8) 



thus defining the coefficients / n ,n' characterizing the dif- 
fusion and providing the S— matrix (see below), the fac- 



3 



tor 1 jr in the second term stemming from the expression 
of a spherical wave in 3D. In a similar way, one can in- 
terpret the transformation Eq. 0} in terms of such a 
S— matrix, and hence to check unitarity. 

In the cosmological one dimensional case, the 1/r fac- 
tor of Eq. (JSJ is absent, the role of the radial variable r 
is played by the absolute value \rj\ of the conformal time, 
and the unit vectors n and n' can point towards only 
two directions, namely +fj and —77, so that n • n' = ±1. 
In what follows, we shall denote the time directions ±77 
simply by the symbols + or — . 

Let us consider the case of the figure on the right, 
i.e. the cosmological case with the potential V u {rj). In 
this case, the analogous form of for u in Eq. © is 
u n n '- After the bounce, the solution for u is 



^aftcr 



1= ( A ni ( 



ik\rj\ 



B„ 



Jk\r,\^ > 



(9) 



from which, in analogy with Eq. |JSj, one gets n = —fj, 
n' = +fj and therefore / |_ is simply the factor B UI / \/2k 



evaluated for A ln = V2fc and B 1 = 0, i.e. / , 

see Eq. Q. Before the bounce, u is given by 



Y*/X, 



A, 



^before 



ife|j7| 



1 



A 



+ (A 1 -A m )e 



ik\r)\ 



(10) 

from which, since one still has n = — fj but n' = —77, 
the comparison with Eq. © provides the coefficient /__, 
namely the value of [A l — A lu ) / \pl~k evaluated for A lu = 
V2k, i.e. /__ = 1/X — l. Reasonning along the same line 
but in the reverse direction for n gives the coefficients 
/++ =1/X-1 and /+_ = -Y/X. 

The total wavefunction \&(M) is the superposition of 
all the solutions for all incident directions weighted by 
arbitrary coefficients F n , i.e. 



F vfr ,rPn 



cx e- ikr F_ a ,-e iKr S n>Il ,F a ,, 



ikr t 



(ii) 



with which one defines the S— matrix as the coefficients 
in front of the outgoing modes. Similarly, in the bouncing 
case, one has the discrete sum 



n— ±r) 



(12) 



leading to 

u+ = F_c-* k M + [f_ + F_ + (1 + /++) F+] e ifc l"l, (13) 

and a similar formula by changing fj —> —fj. Gathering 
all the above results, this leads a two-by-two S— matrix 
that reads 



S = 1 + / = l(>.f), 



(14) 



which, because \X\ 2 — \Y\ 2 — 1 as discussed in the pre- 
vious sections, is unitary (SS^ = 1). We have there- 
fore built, from the non- unitary transformation T u , the 
relevant unitary S— matrix, emphasizing once more the 
probability conservation. 



IV. TIME INVERSION 

In Refs. 0,01, we expanded the matrix T u of Eq. Q in 
the parameter T, which gives a measure of the deviation 
to the null energy condition p+p > at the bounce. For 
T -c 1, the potential V u (rf) goes to a function akin to 
a Dirac 5 distribution, namely V u {rj) — > —CxS(rj), with 
Cf cx Y -1 / 2 3> 1. With such a potential, the perturba- 
tion u crosses the bounce unaltered, i.e. [it] =0, while its 
time derivative jumps by the amount [u'] — — Cxw(0). 
These relations allow to obtain the matrix exactly, and 
we found 



T u = 



/1 Cr 



.C r \ 
~ l 2k * 



1 .Cy 
1 + 1 — 
2k 



= l-i- 



C r ( 1 1 



2k 



-1 -1 



(15) 



which is indeed of the general form obtained in Eq. . 
We use this opportunity to correct a misprint in Eq. (77) 
of Ref. in which the factor i appears in the last term in 
the denominator whereas it should be at the numerator. 

In the presentation made in Ref. 0, we explicitly as- 
sumed the limit n(n + 2)T <C 1 and thus neglected the 
unit matrix part of T u . Using Eq. (|15f) . this gives 



T 



(approx) 



ll 2k 



1 1 

-1 -1 



(16) 



in perfect agreement with Eq. (70) of Ref. Q . As a con- 
sequence, the determinant of the above matrix vanishes 



det 



y(appr 



= 0, 



(17) 



whereas, in fact, the determinant of the exact matrix is 
unity, det (T u ) = 1. 

With a time-symmetric potential V u (—rj) = V(rj) and 
only a second order time derivative term, one expects T u 
to be a time reversal invariant transformation. In other 
words, if the initial conditions are no longer fixed in the 
region I but in the region III, the final result, now ex- 
pressed in region III and no longer in region I, should 
obviously be the same, in particular should still be pro- 
portional to Y -1 / 2 . If one wishes to check this property 
using r p( a PP TOX >^ one f aces the problem that it is not in- 
vertible. Hence, neglecting the term 1 in Eq. I|15fl has led 
to an apparent non symmetrical result while the problem 
clearly possesses this symmetry. The resolution of this 
paradox is obvious: the symmetry has been lost when 
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the approximation has been performed and T u , which 
is invertible since det(T u ) ^ 0, should be used in this 
case. We are only authorized to neglect the term 1 in 
Eq. I|15l) if we want to obtain a close form of the tran- 
sition matrix, as done in Ref. [5j, but not if we want to 
check the time reversal invariance of the matrix. The fact 



that det 



-,(approx) 



= is actually nothing but a consis- 
tency check that the calculation performed in Ref. [B| is 
correct. Indeed, without this property, the inversion of 
^approx) wou y nave b een possible, leading to an effect 
cx T 1 / 2 whereas the time reversal invariance implies that 
it should scale as T -1 / 2 in both time directions. The fact 



that det 



T'V 
u 



(approx) 



can therefore be interpreted as 



the "signal sent by the theory" that we are not allowed 
to use j^ a PP lox ) to check the symmetry of the system. 

It remains to be checked explicitly that the matrix T u 
has the correct time reversal symmetry. For this pur- 
pose, it is sufficient to calculate the inverse of T u and one 
obtains 



T" 1 (C T ) = T u (-C T ) 



(18) 



Therefore, the amplitude of the transition matrix remains 
proportional to T -1 / 2 as required. The only difference 
is that the matrix of the inverse process should now be 
calculated with the parameter — Cx rather than +Cx- 
The origin of this minus sign can easily be understood. 
Indeed, given the perturbation equation time inver- 
sion( in the case at hand, given the analogy with the time- 
independent Schrodinger equation, this operator is simi- 
lar to parity) T manifests itself only in the matching con- 
dition on the first time derivative v! (which is not invari- 
ant) adding an extra minus sign, namely T([it']) = —[it]. 
This is exactly equivalent to changing the sign of Cx- 
Comparing with Eq. I|15|l . we see that 



T[T u (C T )]=T u (-C T ) = T- 1 (Cr). 



(19) 



which merely expresses the familiar result that time in- 
version must satisfy T 2 = T. Note that for an opera- 
tor to be an arbitrary function of some parameter such 
as Cx, the only way to ensure that T 2 = T is that 
T(Cx) = — Cx as indeed is the case here. 

Finally, the explicit form of the S— matrix in the limit 
of Refs. |El3 is 



S = 



1 



1 



iC r 
2k 



( 



X 




-1 

1 o 



(20) 



where in the last approximation we have taken advantage 
of the condition n(n + 2)T < 1 as before. 



V. CONCLUSIONS 

We have presented a brief note on the propagation of 
primordial perturbations through a bouncing phase in 
the early universe to interpret this evolution in terms 
of either transmission and reflexion coefficients across 
the potential of a non-relativistic Schrodinger-like equa- 
tion or in terms of the corresponding S— matrix that 
can be constructed out of the model. We have demon- 
strated that the results obtained in Refs. 0,0 imply that 
R B + T B = 1 and that the S— matrix is unitary, thereby 
checking the mathematical consistency of those results. 
Finally, we have also emphasized that the structure of 
the evolution equation is invariant under time reversal. 
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